Connection between matrix-product states and superposition of Bernoulli shock 

measures 
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We consider a generalized coagulation-decoagulation system on a one-dimensional discrete lattice 
with reflecting boundaries. It is known that a Bernoulli shock measure with two shock fronts might 
have a simple random-walk dynamics, provided that some constraints on the microscopic reaction 
rates of this system are fulflUed. Under these constraints the steady-state of the system can be 
written as a linear superposition of such shock measures. We show that the coefBcients of this 
expansion can be calculated using the finite-dimensional representation of the quadratic algebra of 
the system obtained from a matrix-product approach. 
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I. INTRODUCTION 

It is known that the time evolution of a product shock 
measure in some of the one-dimensional driven-difFusive 
systems is similar to that of a random- walker, provided 
that the microscopic reaction rates of the system lie on 
a certain manifold These systems can be defined 

both on infinite lattices or finite lattices with boundaries. 
For each of these systems one can define a microscopic 
shock position and calculate the exact hopping rates of 
the traveling wave in terms of the transition rates of the 
microscopic model. The existence of such processes im- 
plies a rather remarkable property. Shocks behave like 
collective single-particle excitations on the lattice scale 
which results in reduction of the exponentially large num- 
ber of microscopic internal degrees of freedom to an only 
polynomially large number of macroscopically relevant 
degrees of freedom. Since the shock position moves like 
a biased single-particle random-walk, the shock measure 
evolves in time into a linear combination of shock mea- 
sures; therefore, a linear combination of shock measures 
may be a stationary measure. 

The stationary measure of some of the one-dimensional 
driven-difFusive systems can also be obtained using a 
matrix-product approach (for a review see 0]). Accord- 
ing to this approach the stationary probability of a given 
configuration is written in terms of the expectation value 
of a product of non-commuting operators associated with 
different states of each lattice site. In the case of nearest- 
neighbor interactions these operators should satisfy a 
quadratic algebra. The matrix representations of some 
of these quadratic algebras have already been studied in 
different of cases (see Q and references therein) . 

Since the steady-states of these systems are unique, 
one can ask about the relation between these two ap- 
proaches. It has been shown that if the steady-state of a 
system with nearest-neighbor interactions can be written 



in terms of superposition of product shock measures with 
a single shock front, then a two-dimensional representa- 
tion of the quadratic algebra would be enough to write 
the steady-state of the system as a matrix-product state 
Q. In this case the conditions, under which a product 
shock measure with a single shock front has a random- 
walk dynamics, are exactly those for the existence of a 
two-dimensional matrix representation of the quadratic 
algebra. A couple of examples have already been studied 

in a. 

In a recent paper Q the authors have studied a gener- 
alized non-conservative reaction-diffusion system defined 
on a finite lattice with reflecting boundaries with the fol- 
lowing non-zero reaction rates: 



A 
A 
A 

0H 
A- 



- A ^ A + (d with rate UJ32 
h ^ + A with rate ^23 
h A ^ A + with rate 0)34 
h A ^ + A with rate UJ24 

- A ^ A + A with rate UJ42 
h ^ A + A with rate UJ43 



(1) 



'Electronic address: ' farhad@ipm.ir| 



in which A and stand for the presence of a particle 
and a hole in each lattice site respectively. Note that 
there is no injection or extraction from the boundaries. 
It has been shown that the steady-state of this system 
can be written as a matrix-product state under some con- 
straints @. The quadratic algebra of the system has a 
four-dimensional representation in this case. The time 
evolution of a product shock measure with two shock 
fronts has also been studied for this model. It has been 
shown that the shock positions in this measure have sim- 
ple random-walk dynamics provided that the same con- 
ditions for the existence of the four-dimensional matrix 
representation for the quadratic algebra of the system are 
fulfilled by the microscopic reaction rates 0]. In other 
words these shock distributions form an invariant sector 
under the time evolution of the system i.e. a shock mea- 
sure evolves into a linear combination of shock measures 
with different shock positions. In this case one can write 
the steady-state of the system as a linear superposition of 
shocks; however, this has not been shown yet. The com- 
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plete phase diagram of this system has also been studied. 
It has been found that it has two different phases depend- 
ing on the hopping rates of the shock positions. 

In this paper we aim to study the steady-state of the 
exactly solvable system defined in ([T]) using a different 
approach than the matrix-product approach. Wc will 
show that the steady-state of the system can be ex- 
pand in terms of product shock measures with two shock 
fronts. We will specifically show that the coefficients of 
the expansion can be uniquely calculated using the four- 
dimensional matrix representation of the quadratic alge- 
bra introduced in In the following we will start with 
the mathematical preliminaries and bring the time evo- 
lution equations for a product shock measure with two 
shock fronts generated by the Hamiltonian of the system. 
By considering a linear superposition of such shock mea- 
sures we construct the steady-state of the system and 
then calculate the coefficients of this expansion. We will 
provide the summery and discussion at the end of the 
paper. 



with reflecting boundaries the Hamiltonian H is of the 
following form: 



L-l 

H — ^ hkM+i 



(3) 



fe=i 



in which: 



where T is a 2 x 2 identity matrix and /i is a 4 x 4 
matrix representing the bulk interactions. In the basis 
(00, 01, 10, 11) with the following basis vectors: 



= 1) = ( J ) , = 0) = J ) 



for /c = 1 , • • • , L 
(4) 

the Hamiltonian h for the system defined by H]) can be 
written as: 



II. TEMPORAL EVOLUTION OF SHOCKS 

The microscopic state of a one-dimensional reaction- 
diffusion system may be described by a set of occupation 
numbers {ni, ■ • • ,71^} where = 0, 1 is the number of 
particles on site k on a lattice of L sites. The time evolu- 
tion of the distribution P({ni, ■ ■ ■ , n^}; i) in this system 
is defined by a continuous-time master equation which 
can be written in terms of the quantum Hamiltonian for- 
malism. The stochastic Hamiltonian H, whose matrix 
elements are the transition rates between configurations, 
generates the time evolution. The Markovian time evo- 
lution can be written in the form of an imaginary time 
Schrodinger equation 



-\P(t))=H\Pit)). 



(2) 



For a single-species system with nearest-neighbors inter- 
actions defined on a one-dimensional lattice of length L 
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(5) 

In [9] the authors have shown that a product shock mea- 
sure with two shock fronts defined as: 



1 X ^n — m—l / 1 \ ®L — n+1 

^ 



(6) 

for < TO < n — 1 and 1 < n < L + 1 might evolve 
according to 2-particle random-walk dynamics provided 
that: 



1 — P ^2i + ^34 _ ^^23 _ ^^32 

p UJ42 + ^43 ^43 ^42 



(7) 



Note that in ([6]) two auxiliary sites and L+1 are defined 
for convenience. The time evolution of the product shock 
measure ^ is obtained to be: 



H\Pm,n) = Sir\Prn+l,n) + <5l/ |P„i-l,„) + S2r\Pm,n+l) + S2l\Pm,n-l) ~ (<5lr + <5l/ + i2r + ^2;)l^m,n) 

for TO = 1, • • • , L — 2 and n = m + 2, • • • , L 

H\Po^n) - -S\Pi,n) + S2r\Po,a+l) + ^2;|^0,«-l) " {^6 + 52r + <52/)|Po,n) for n = 2, • • • , L 

i^|Pm,L+i>='5i,|P,„+i,L+i>+5i,|P„_i,L+i)+5|P„,,i)-(,5i, + ,5i,+5)|P„,.i+i) for TO = 1,---,P-1 

H\Po,L + l) = -S\Pl.L + l) + S\Po.l) 

H\Pm,ni+i) =0 for m = 0, • ■ • , L 

I 

in which we have defined: As can be seen in the bulk of the lattice the shock posi- 

Su ■= — tions move like two biased lattice random- walks; there- 

dir := UJ43 p + ^24P, 

S2I := •^42'^—^ + ^34/3, (9) 
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fore, provided that we are far from the boundaries, the 
velocity, and also the diffusion coefficient of each shock 
front, can be easily calculated from ([8]). The hopping 
rates of the left (right) shock front are Sir,Sii {S2r,S2i)- 
The velocity of an individual shock front is then Vi — 
Sir — Sil for i — 1,2 and also the diffusion coefficient is 
A-(^^r + <5^^)/2for^-l,2. 



III. STEADY-STATE OF THE SYSTEM 

The last equation in ([8]) implies that an empty lattice 
is a trivial steady-state for the system; however, one can 
construct a non-trivial steady-state by considering a lin- 
ear superposition of |P„i,ti}'s as follows: 



1 



L L+1 



m— n— m-t-1 



'4^m,n I Pm^n) 



and find the coefficients ipm,n's by requiring that: 



HIP* 



0. 



(10) 



(11) 



The probability of finding the system in the configuration 
{n-i, ■ ■ ■ ,71^} is now given by: 



P*{{nu--- ,nL}) = {{ni 



{nL\)\P*) (12) 



in which \nk) ior k — 1, ■ ■ ■ , L is defined in Note that 
there are L+1 states \Pk^k+i)^s ior k ~ 0, ■ ■ ■ , L which 
all point to an empty lattice and should be considered as 
a single state. The coefficient of this state in ([TOl) will be 
called ip' . 

It is a lengthy exercise, but nevertheless straight- 
forward to find the equations governing V'm.n's using ([S]), 
(fTU)l and It turns out that these coefficients should 
satisfy a set of fifteen equations which are given in Ap- 
pendix. In order to find the non-trivial steady-state of 
the system one should solve these equations and calcu- 
late the coefficients V'm.n's. These equations have already 
been solved in [4| for the following special tuning of the 
parameters: 



(13) 



In this paper we are not going to solve these equations 
directly, but instead, we will find their solutions using 
the results obtained from the matrix-product approach. 

The non-trivial steady-state of the system has already 
been obtained using the matrix-product approach [9|. 
Let us assign the two operators D and E to the exis- 
tence of a particle and a vacancy at each lattice site re- 
spectively. According to this approach the normalized 



matrix-product steady-state is given by: 

\pi^^^{{w\{iy''^ \v)) (14) 

in which Zl = {{W\{D + E)^\V)). In contrast to (§], 
here we use a new notation for the two vectors \V)) 
and {{W\ to emphasis that the vectors in the configu- 
ration space denoted by | • • • ) are different from these 
two vectors which exist in an auxiliary space. Using (jl2p 
the probability of finding the system in the configuration 
{rii, • • ■ tUl} is now given by: 

P*({ni, • • • ,nL]) = — ((M/| WinuD + (1 - nk)E)\V)). 



k=l 



In the authors have found that the two operators D 
and E besides the two vectors \ V)) and {(W\ have a four- 
dimensional matrix representation provided that the con- 
straints ([7]) are satisfied. Because of uniqueness of the 
steady-state of the system, the expression pUj) should be 
equal to the expression p4|) . This will help us calculate 
the coefficients V'm.n's. It can easily be checked that the 
results obtained in this way satisfy the equations govern- 
ing V'm.n's. 



IV. EQUIVALENCE OF TWO DIFFERENT 
APPROACHES 



Let us define a new measure: 



1 

-P 



\Pm,n) — I 1 — p 

/ \ P / \ -o / 

(15) 

for < TO < n — 1 and 1 < n < L + 1 which is orthog- 
onal to the product shock measure ([6]) according to the 
following rule: 

{Pyn'^n' \Pm^n} Syyi ^n'Sn^n'- (1^) 

Using this new measure and (fTU]) one finds: 



'^m^n — ^L{Pm,n\P 



(17) 



for < TO < n — 1 and 1 < n < L + 1 which uniquely 
determines the coefficients '4'm,nS- Since the two steady- 
states obtained from the superposition of the product 
measures and the one obtained from the matrix-product 
approach should be equal, one finds from and 
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„ = {{W\{E - l^D)"'{^DY'-"'-^{E - ^D)^-^'+^\V)) for < m < L - 1 and m + 2 < n < L + 1, 

(18) 

^'^{{W\{E-'-^Df\V)). 



If there exists a matrix representation for the quadratic 
algebra of the system, one can calculate the coefficients 
V'm.n's using (jlSp . and it means that the matrix-product 
steady-state of the system can actually be expand in 
terms of a linear superposition of product shock mea- 
sures with two shock fronts. 



Using the diagonal four-dimensional matrix represen- 
tation of the quadratic algebra of the system, first intro- 
duced in 9'], we have calculated the coefficients ipm.n's. It 
turns out that these unnormalized coefficients are given 
by the following expression: 



V'. 



f .^lr(£li+£2T0_-,5„ 



( a2((-5i;+*2r) 



— (t^)'"(#^)"(( x""^x" - (1^41^)) for < m < L - 1 and m -I- 2 < n < L + 1, 



.1,1 _ pSuSir I S2r(S2l+&lT) \2 I ('?21+'?lr)p ( hjL ( ^ — n)'^'\^ 1 (ilrS2r\L 1 

''^ ~ {Su+S2,-)^SuS2l-Sl^S2j ^{Sl,.-5u(l-py')iS2l-S2Al-py')^S2r f'' (Slr-Sud-P)') ^ SuS2l ' (52i-52.(l-p)^) 



(19) 



in which p = \ — {5ir + 621) / 52r)- These coefficients 
are functions of the four parameters Ji,-, 5ii, and S21 
as one should expect. As we mentioned it can be easily 
verified that 4'm,nS besides '0' in satisfy the equa- 
tions (jAl[) . Note that in order to calculate the coefficient 
-0' in one should also include the fact that an empty 
lattice (the trivial steady-state) should be excluded by 
applying the constraint P*({0, • • • , 0}) oc {{W\E^\V)) = 
0. One can also check that for the special tuning of the 
parameters given in (jl3p . the unnormalized coefficients 
(fT9| reduce to the ones calculated in [4] (up to a multi- 
plicative constant). 

As we mentioned the system has two different phases 
which are separated by a coexistence line. On this coex- 
istence line, where one of the shock fronts performs an 
unbiased random-walk, the matrix representation of the 
quadratic algebra cannot be diagonalized and one should 
use a non-diagonal matrix representation. It can be eas- 
ily shown that even in this case one can calculated V'm,n's 
and ■(/'' using (|18p and the non-diagonal matrix represen- 
tation of the quadratic algebra introduced in Q . 

V. CONCLUSION 

In this paper we have studied the steady-state of a gen- 
eralized coagulation-decoagulation system defined on a fi- 
nite lattice with reflecting boundaries. The steady-state 
of this system has already been found using a matrix- 
product approach; however, in this paper we have shown 
that the steady-state of the system can equivalently be 



constructed by considering a linear superposition of prod- 
uct shock measures with two shock fronts. The key point 
is that the conditions under which the shock fronts have 
simple random- walk dynamics, are exactly those neces- 
sary for the existence of a four-dimensional representa- 
tion for the quadratic algebra of the system. This is 
quite non-trivial since it has not been generally shown 
that the existence of a finite-dimensional matrix repre- 
sentation for the quadratic algebra of a given system is 
an indication that the steady-state of the system can be 
written as a superposition of Bernoulli shock measures. 
The results of the study of this exactly solvable system 
besides those obtained in [1] might bring us to this con- 
clusion: that the existence of a finite-dimensional ma- 
trix representation for the algebra of a one-dimensional 
reaction-diffusion system is a signal that the steady-state 
of the system can be written as a superposition of prod- 
uct shock measures; nevertheless, the proof still remains 
as an open problem. 



APPENDIX A: EQUATIONS GOVERNING Vm n'S 

As we mentioned, the non-trivial steady-state of the 
system can be written as a linear superposition of 
Bernoulli shock measures with two shock fronts of type 
By requiring that PH)) is the steady-state of the 
system and using ([5]), pU|) and pTjl one finds that the 
coefficients of this expansion should satisfy the following 
difference equations: 
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{S_ - bll)^Q.2 - {8ir + ?>2l) Sm=l '0m,m+2 - ((^ + (5lr)^L-l,L+l = 
((5 - ((52r + '520)V'O,2 + '52i'/'0,3 = 

^2rf/'0,n-l + ((5 - {hr + ^zOV'O^n + ^2iV'0,n+l + (^IfV^l,™ =0 for n 2, L 
<5?/'0,L + l + hri'a,L-\ + - (<^2r + ^2l))i'a.L + (^iz-^i^L = 
^2r'0O,L + <5ijV'1,L+1 = 

^V'O.S + (l5lr- + (^IZ + hr + '52;)V'1.3 " hl^^XA = 

5V'0,n - 52rV'l,n-l + ((^Ir + + fe- + (52/)'(/'l,n " <52;l/'l,«+l " Sll1p2,n = for n 7^ 3, X 
^(■00, L - V'l^L+l) - S2ri>l,L-l + [Sir + Su + S2r + S2l)i>l,L " Sil'lp2,L = 
Slpo.L + l + {S + Sir + du)'lpi,L + l - 5u1p2.,L + l - '52rV'l,L = 

Slr^^^Jn-l,L + l - {5 + Sir + Sll)i>m,L + l + Sui'm+l.L + l + S2r1prn,L = foi TTl ^ 1,L - 1 
5_lri^L~2X + l - [5 + Sir + (5l/)'(/'L-l,L + l = 

S_^m,L + l + Slr1pm-1,L + S2r'4^mX-l - {Sir + Su + S2r + S2l)lpm,L + Snlpm+l.L =0 for m ^ 0, 1 , L - 1 , L - 2 
SlpL-2,L + l + Sir-ipL-dX - (Sir + Su + S2r + S2l)ipL-2X = 

Sir1pm-l,n - [Sir + Su + S2r + 52l)ll^m,n + ^2;l/')n,n+l =0 for m 7^ 0, 1, L - 2, L - 1 and n = TO + 2. 
Slri'm-l,n + <52r'0m,n-l " {Sir + Su + S2r + S2l)tpm,n + S2l1pm,n+1 + '^uV'm+l.n = 

for TO 7^ 0, 1, L — 1 and n^m + 2,L,L + l 

I 

Note that S is related to the shock hopping rates through 
the last relation in 
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